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Abstract— Based on the Ogorodnikov-Milne model, we analyze the proper motions of 95 
633 red giant clump (RGC) stars from the Tycho-2 Catalogue. The following Oort constants 
have been found: A = 15.9 ± 0.2 km s~^ kpc~^ and B = —12.0 ± 0.2 km s~^ kpc~^. Using 
3632 RGC stars with known proper motions, radial velocities, and photometric distances, 
we show that, apart from the star centroid velocity components relative to the Sun, only the 
model parameters that describe the stellar motions in the XY plane differ significantly from 
zero. We have studied the contraction (a negative K effect) of the system of RGC stars as a 
function of their heliocentric distance and elevation above the Galactic plane. For a sample 
of distant (500-1000 pc) RGC stars located near the Galactic plane (|2'| < 200 pc) with an 
average distance of d = 0.7 kpc, the contraction velocity is shown to be Kd = —3.5 ± 0.9 
km s~^; a noticeable vertex deviation, l^y = 9° ± 0.5°, is also observed for them. For stars 
located well above the Galactic plane {\Z\ > 200 pc), these effects are less pronounced, 
Kd = —1.7 ± 0.5 km s~^ and l^y = 4°. 9 ± 0°.6. Using RGC stars, we have found a rotation 
around the Galactic X axis directed toward the Galactic center with an angular velocity of 
—2.5 ± 0.3 km s^^ kpc^^, which we associate with the warp of the Galactic stellar-gaseous 
disk. 

INTRODUCTION 

The red giant clump (RGC) stars are distributed fairly uniformly over the celestial sphere. 
Since they occupy a compact region on the Hertzsprung-Russell diagram and, hence, have the 
same luminosity, reliable photometric distance estimates can be obtained for them. These 
stars are of great interest in analyzing the three-dimensional kinematics of various Galactic 
subsystems. 

First of all, we are interested in such peculiarities as the expansion/contraction of stellar 
groups, the deviation of their vertex from the direction of the kinematic center, i.e., the 
Galactic center, and the connection of the stellar kinematics with the warp of the Galactic 
stellar-gaseous disk. 

Various authors have found a negative value of the Oort parameter K responsible for the 
expansion/contacrtion. This means that the system of stars is in the state of contraction. 
We are talking about the plane effect that describes the contraction in the Galactic XY 
plane. The negative K effect in the motion of distant OB stars was detected by Torra et al. 
(2000) and Bobylev (2004). It was also revealed in the motion of open star clusters (Bobylev 
et al. 2007) and in the motion of dwarf stars (Rybka 2004a) and giant stars (Rybka 2004b, 
2008). 



At present, the nature of the negative K effect is not completely clear. It is hypothesized 
that this effect is related to: (i) the pecuharities of stellar radial velocity measurements 
(Pont et al. 1994), (ii) the influence of the bar at the Galactic center (Point et al. 1994), 
and (iii) the influence of the Galactic spiral structure (Rohlfs 1977; Metzger et al. 1998). 
In particular, Fernandez et al. (2001) and Bobylev et al. (2006a) made an attempt to take 
into account the influence of a spiral density wave; it failed, because the negative K effect 
was retained. 

Rybka (2008) showed that the RGG stars have different vertex deviations and differ 
significantly by the magnitude of the K effect, depending on the Galactic latitude. In this 
respect, it is interesting to trace the changes in kinematic parameters of these stars located 
in various layers parallel to the Galactic plane. 

Analysis of the large-scale structure of neutral hydrogen has shown that the gaseous 
disk in the Galaxy is warped (for a review, see Burton, 1988). The results of studying this 
structure on the basis of currently available data on the HI distribution arc presented in 
Kalberla and Dedes (2008). This structure is revealed by the spatial distribution of stars 
and dust (Drimmel and Spergel 2001) and by the distribution of pulsars in the Galaxy 
(Yusifov 2004). Having analyzed nearby Hippacros (1997) stars, Dehnen (1998) showed that 
the distribution of their residual velocities Vy — Vz was in satisfactory agreement with various 
rotation models of the warped disk. Miyamoto et al. (1993) and Miyamoto and Zhu (1998) 
determined the rotation parameters of the warped stellar-gaseous disk by analyzing giant 
stars of various spectral types and samples of Hipparcos 0-B5 stars. Thus, there is positive 
experience in solving this problem using data on stars relatively close to the Sun. 

The goal of this paper is to study the kinematic peculiarities of a wide solar neighborhood 
on the basis of data on the RGG stars from the Tycho-2 Gatalogue (Hog et al. 2000) selected 
by Gontcharov (2008) and, in particular, to establish the connection of the kinematics of 
RGG stars with the warp of the Galactic stellar-gaseous disk. 

1 DATA 

At present, two lists of candidates for RGG stars compiled from a combination of Tycho-2 
and 2MASS data (Skrutskie et al. 2006) have been published. These are the catalogs by 
Rybka (2006, 2007) and Gontcharov (2008). 

The catalog by Rybka (2007) contains about 60000 RGG candidates from Tycho-2 
with photometric characteristics from 2MASS with distances > 100 pc. The catalog by 
Gontcharov (2008) contains about 97000 stars with data from the same catalogs as those of 
Rybka but without any distance constraint. Therefore, the number of stars is larger. In con- 
trast to the catalog by Rybka (2007), the catalog by Gontcharov (2008) gives a photometric 
distance estimate for each of the stars corrected for interstellar extinction. In this respect, 
the catalog by Gontcharov (2008) is more convenient for solving the kinematic problems 
considered here. 

The radial velocities taken from the PGRV (Gontcharov 2006; Bobylev et al. 2006b) and 
RAVE (Steinmetz et al. 2006) catalogs are given for 4163 stars. These are relatively close 
stars with known Hipparcos trigonometric parallaxes. 

Figure 1 gives an idea of the three-dimensional spatial distribution of stars. The left, 
middle, and right columns present the distributions of stars in the XY, XZ, and YZ planes, 
respectively. The rows from top to bottom provide the distributions for the complete sample 



of stars and the samples of stars with distances d < 400 pc, 400 < d < 600 pc, and 
600 < d < 1600 pc. The distance ranges were chosen in such a way that each sample 
contained approximately the same number of stars. The figures reveal some structural 
peculiarities of the samples. Thus, for example, clumps inclined by ~ 20° to the X axis, 
typical of young stars associated with the Gould Belt, are clearly traceable in the ZX plane 
(c) and (d). 



In this paper, we use a rectangular Galactic coordinate system with the axes directed away 
from the observer toward the Galactic center (/ = 0°,b = 0°, the X axis or axis 1), along 
the Galactic rotation (/ = 90°, b = 0°, the Y axis or axis 2), and toward the North Galactic 
Pole (6 = 90°, the Z axis or axis 3). 

2.1 The Ogorodnikov-Milne Model 

Using the three-dimensional velocity field. In the linear Ogorodnikov-Milne model, 
we adhere to the notation that was used by Clube (1972, 1973). The observed velocity V(r) 
of a star with a hehocentric radius vector r is described, to the terms of the first order of 
smallness r/i?o <^ 1, by the equation in vector form 



where 'Vq{Xq, YqjZq) is the velocity of the Sun relative to the centroid of the stars under 

consideration and V is the residual velocity of the star. Here, the residual stellar velocities 
are assumed to have a random distribution. M is the displacement matrix whose components 
are the partial derivatives of the velocity u(mi, U2, M3) with respect to the distance r(ri, r2, r^), 
where u = V(i?) — V(i?o), while R and Rq are the Galactocentric distances of the star and 
the Sun, respectively. Then, 



All nine components of the displacement matrix can be determined if the distances, radial 
velocities, and proper motions of the sample stars are known. In this case, the conditional 
equations are 



2 THE METHODS OF ANALYSIS 



V(r) = V© + Mr + V, 



(1) 




(2) 



Vr — —Xq cos b COS / — 
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+ cos b sin b cos /M13 + cos^ b sin / cos /M21 + 
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+ sin2 6M33), 
4.7Arfii cos b = Xq sin / — Yq cos / + 
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— sin b sin /M13 + cos b cos^ /M21+ 
+ cos b sin / cos IM22+ 
+ sin 6 cos IM23), 
i.Tirfib = cos / sin 6 + 
+Yq sin Z sin 6 — Zq cos b+ 
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— sin b cos 6 sin I cos ZM12 — 
— sin^ b cos ZMi3 — sin b cos 6 sin I cos ZM21 — 
— sin b cos 6 sin^ ZM22 — sin^ b sin IM23+ 
+ cos^ 6 cos ZM31 + cos^ 6 sin ZM32+ 
+ sin 6 cos 6M33), 



(5) 



where the stellar proper motion components are in mas yr^^ (milliseconds per year), the 
radial velocity Vr is in km s"\ and the heliocentric distance to the star r is in kpc; we 
either calculate the latter based on the known trigonometric parallax r — I/tt or take the 
photometric distance dphot- Equations (3)-(5) contain twelve sought-for unknowns - three 
components of the velocity V0(X0, 1^0,^0) in km s^^ and nine components of Mpq in km 
s~^ kpc^^. The system of equations (3)-(5) is solved by the least squares method. 

The Oort constants can be found as follows: A = 0.5(Mi2 + M21), B = 0.5(M2i - Mu), 
C — 0.5(Mii — M22) and K — 0.5(Mii + M22); we determine the vertex deviation Ij^y from 
the relation tan tan2/a;y ~ —C/A. The angular velocity of Galactic rotation can be found 
from the relations {^z)Ra = B — A = — M12. The equations are written in such a form that 
the negative sign of the angular velocity corresponds to the direction of Galactic rotation. 

The matrix (tensor) M can be divided into symmetric, M"*", and antisymmetric, M~, 
parts. Following Ogorodnikov (1965), we call them the local deformation tensor and the 
local rotation tensor, respectively: 



Using only the proper motions. When only the stellar proper motions are used, 
one of the diagonal terms of the local deformation tensor is known to become uncertain. 
Therefore, we can determine only the differences, for example, (M+ — M^^ and (M3+ — M^^. 
In this case, we use the conditional equations 




2 \ dvq drpj 
(p,g = l,2,3). 



(6) 



A.TArjii cos b = Xq sin / — Yq cos / — 



(7) 



r{M^^ cos / sin 6 — M^g sin / sin 6 + M^^ cos b+ 



+Mj^ cos 21 cos b — sin / sin b+ 



+M+ cos / sin 6 - 0.5(M+ - M+) sin 21 cos b), 
4:.74rfib = Xq cos / sin b + Yq sin / sin b — 



—Zq cos b + r (M32 sin / — cos /— 



-0.5M+ sin 21 sin 2b + M+ cos / cos 26+ 
+M+3 sin / cos 2b - 0.5(Af + - M+) cos^ / sin 26+ 
+0.5(M3+-Af+)sin26). 

is equivalent to the Oort constant B. In accordance with the chosen coordinate system, 
the positive directions of rotation are: from axis 1 to 2, from axis 2 to 3, and from axis 
3 to 1. M+ is equivalent to the Oort constant A. The diagonal components of the local 
deformation tensor M+, M^^, M3+ (by definition, they coincide with Mn, M22, M33) describe 
the general contraction or expansion (depending on the sign) of the entire stellar system. 
Thus, the system of conditional equations (7)-(8) includes eleven sought-for unknowns to be 
determined by the least squares method. 

Although the unknowns and M+ as well as and M^^ in Eq. (7) being determined 
cannot be separated between themselves, because they have the same coefficients (sin I sin b 
and cos I sin b, respectively) , all variables can be separated through the simultaneous solution 
of the system of equations (7) and (8). 

The Oort constant K = (M+ + M^^/2 can be calculated on the basis of the quantities 
found by solving the system of equations (7)-(8) as follows: 

K = ((M^ - M+) - 2(M3+3 - M+))/2 (9) 

by assuming M3+ to be close to zero. Therefore, we first determine M3+ by solving the 
system of equations (3)-(5) based on the sample of RGC stars with known trigonometric 
parallaxes, radial velocities, and proper motions to test this assumption. We determine the 
vertex deviation l^y from the relation tan 2/^,2^ — —CjA — — 0.5(M+ — M^^jM^^. 

2.2 The Statistical Method 

We use the well-known statistical method (Trumpler and Weaver 1953; Parenago 1951, 1954) 
to determine the parameters of the residual velocity (Schwarzschild) ellipsoid. It consists 
in determining the symmetric tensor of moments or the tensor of residual stellar velocity 
dispersions. When simultaneously using the stellar radial velocities and proper motions 
to find the six unknown components of the dispersion tensor, we have six equations for 
each star. The semiaxes of the residual velocity ellipsoid, which we denote by 0-1,2,3 can be 
determined by analyzing the eigenvalues of the dispersion tensor. We denote the directions 
of the principal axes of this ellipsoid by /i,2,3 and 61,2,3- 

3 RESULTS 

The results of solving the system of equations (3)-(5) that we obtained based on the sample 
of 3632 RGC stars with known trigonometric parallaxes, radial velocities, and proper motions 
are presented in Table 1. The kinematic parameters were found using both trigonometric 
parallaxes and photometric distances with the constraints c^/tt < 1 and tt > 1 mas. A 
preliminary analysis shows that we need these constraints to eliminate the stars that spoil 
the statistics and affect significantly the parameters being determined. 

As we see from Tabic 1, M33 is found to be statistically insignificant. Therefore, we have 
every reason to neglect M33 below when determining the K effect from Eq. (9) based only 
on the stellar proper motions. 



The parameters of the Ogorodnikov-Milne model obtained on the basis of 95633 RGC 
stars are presented in Table 2. Based on the results of Gontcharov (2008) and the data of 
Fig. 1, we may conclude that the spatial distribution of very distant RGC stars is highly 
nonuniform. Using RGC stars farther than 1 kpc can lead to undesirable selection effects. 
Moreover, as was shown by Oiling and Dehnen (2003), the nonuniformity of the distribution 
in Galactic longitude affects noticeably the determination of kinematic parameters, in par- 
ticular, the Oort constants. Therefore, we consider stars closer than 1 kpc. Columns nos. 
1-3 in Table 2 give the parameters obtained with various constraints on the heliocentric 
distance of a star d{dphot) columns nos. 4 and 5 give the parameters for stars with distances 
of 500-1000 pc located near the Galactic plane (|Z| < 200 pc) and higher (|Z| > 200 pc), 
respectively. 

We determined the parameters of the residual velocity ellipsoid based on the statistical 
method both for all stars of the catalog and for individual samples. The first column in 
Tabic 3 gives the parameters of the residual velocity ellipsoid for the sample of 3632 stars 
considered above. We estimated the errors in the parameters obtained from this sample 
based on the well-known method by Parenago (1951) using the fourth-order moments. 

The complete sample of stars was divided into individual groups by the physical properties 
of the stars themselves as well as by the heliocentric distance d and the elevation above the 
Galactic plane \Z\. The results are presented in columns nos. 2-4 of Table 3 in the former 
case and in Table 4 and Fig. 2, which shows the residual velocity dispersions as a function 
of |Z| found by the statistical method using the proper motions of 95633 RGC stars, in the 
latter case. 

We divided the RGC stars by their physical properties based on the selection parameters 
"Sell" and "Sel2" listed in the last columns of the catalog by Gontcharov (2008). The 
first parameter points to the displacement of a star from the center of the distribution 
of young RGC stars on the {B — l/)o-reduced proper motion diagram, while the second 
parameter points to the displacement of a star from a similar diagram for old RGC stars 
(Gontcharov 2008). We formed three groups: group (1) consists of stars with a difference 
< AS, where AS = Sell-Sel2, which includes old normal-metallicity stars; the stars of 
group (2) have —0.38 < AS < and belong to old low-metallicity stars; finally, group 
(3) with —1.0 < AS < —0.38 contains young normal metallicity stars. Consequently, the 
third and second groups include the youngest and oldest stars, respectively. The increase 
in average \Z\ as we pass from the relatively young groups to the older ones: \Z\^^^ — 183 

pc, 1-^1(1) = 201 pc, 1^1(2) = 221 pc, indicates that this division identifies the stars by 
age statistically correctly. Since only the proper motions of RGC stars were used in this 
case (their space velocities U,V,W cannot be determined), we estimated the errors in the 
parameters of the residual velocity ellipsoid given in column no. 5 of Table 3 statistically, 
based on the Monte Carlo method. 

Figure 3 shows the distributions of residual UV velocities smoothed by the method de- 
scribed by Bobylev ct al. (2006b) for groups (1), (2), and (3) (separated by the parameter 
?S); we used only 3632 RGC stars with known trigonometric parallaxes, radial velocities, 
and proper motions (see Table 1). As we see from the figures, in comparison with sample 
(3) of "young" stars, the UV distributions for samples (1) and (2) of "old" stars exhibit 
a larger velocity dispersion manifested as the "spreading" of peaks and the appearance of 
new clumps, which is caused by the prolonged action of some dynamical factor (a bar, a 
spiral density wave, and the like). However, it should be noted that the velocity dispersions 
calculated using the known trigonometric parallaxes rather than the photometric distances 



(as with the complete sample) show the highest values for sample (1), not for (2), which is 
in confiict with the results of Table. 3. Quahtatively, the older age of the stars from group 
(1) follows from a comparison of the C/y- velocity distributions shown in Fig. 3. 

Consequently, the method of selection by the parameter AS* we used is a reliable method 
for identifying only the comparatively "young" fraction of stars. A reliable separation of stars 
with various ages among the remaining "old" fraction of stars is unlikely to be possible, as 
suggested by the insignificant difference in velocity dispersions for groups (1) and (2), while 
the velocity dispersions for the young fraction of stars differ significantly. In addition, we 
obtained contradictory results of the selection by the parameter AS for the entire sample 
of ^97 000 stars and the sample of 3632 stars with known parallaxes, radial velocities, and 
proper motions. Clearly, groups (1) and (2) are actually a mixture of stars with various ages. 

In Fig. 4, the contraction velocity Kd is plotted against the heliocentric distance d. 

The parameters being determined depend on the division of the original sample into 
subsamples according to the measured distances. Since all distances contain errors, on 
average, they show systematic deviations from the true distances. These deviations correlate 
with the measured distances themselves (the largest measured distances are, on average, 
overestimated, while the smallest ones are, on average, underestimated). 

Therefore, to ascertain the dependence of the kinematic parameters being determined 
on the distances including the measurement errors, wc performed statistical Monte Carlo 
simulations. To this end, we solved the system of equations (7) and (8) with artificially 
introduced random errors in the distances distributed normally. The division into subsamples 
was made by the derived model distances. We analyzed all stars of the catalog without any 
constraints, but the star was rejected in the case where a negative distance emerged. The 
equations were solved for 100 realizations for each subsample. Wc considered three cases: (i) 
without introducing any errors, so that these results are a good addition to the data of Table 
2; (ii) the model errors are 20% of the measured distance, which is a typical error in the 
photometric distance; (iii) the model errors are 60% of the measured distance (considered 
as the extreme case). 

The simulation results are presented in Figs. 5, 6, and 7, where all eleven parameters of 
the Ogorodnikov-Milne model are plotted against the distance. Figures 5a, 6a, and 7a show 
the velocity components of the Sun relative to the centroid of the stars under consideration 
Yq,Xq; Figs. 5b, 6b, and 7b show the parameters Mpq. 

Analysis of these data leads us to conclude that: (a) when stars with large distances are 
used, a significant bias is observed only for the velocity of the Sun Y©, while all three velocity 
components of the Sun are subjected to a similar bias for a large range of specified errors 
(Fig. 7a); (b) the parameters Mpg being determined are smoothed out as the level of model 
errors increases, which is particularly clearly seen from the middle part of Fig. 7b; (c) it 
is necessary to reject very close {d < 200 — 250 pc) stars for a reliable determination of all 
model parameters, except the velocity components of the Sun, which, on the contrary, are 
determined more reliably from the samples of nearby stars; (d) three parameters are least 
subjected to any biases: M+, M^; (the Oort constants), and (M+ — M^^). 

We calculated the correlation coefficients (coefficient k) between the parameters being 
determined, which arise because the system of equations (7)-(8) is ill-conditioned: (i) k — 
0.44 for the pair {M+) - {M^, (ii) k = 0.41 for the pair {M+ - M+) - (Af+ - Af+), (iii) 
for the sample of most distant stars (Figs. 5-7), k = 0.37 for the pair (M+) — {M^J, while 
for closer stars, the correlation coefficient k is nearly zero for this pair. The correlations 
between the velocity Yq, which is subjected to the largest bias, and all the remaining model 



parameters are nearly zero. 

The slight increase in velocity Yq with distance that we see in Fig. 5a can be explained 
by an asymmetric drift (Bobylev and Bajkova 2007). Indeed, the mean value of AS changes 
gradually from 0.1 (d = 150 pc) to -0.1 (d = 900 pc). 

4 DISCUSSION 

Galactic Rotation Parameters As can be seen from the first column of Table 2, the Oort 
constants are: A — 15.9 ± 0.2 km s~^ kpc~^ and B — —12.0 ± 0.2 km s~^ kpc~^. Based on 
these parameters, we estimate the circular velocity of rotation of the solar neighborhood to 
be Vo = \Ro^z\ = 222 ± 2 km for Rq — 8.0 kpc and the period of its revolution around 
the Galactic center to be 221 Myr. 

The Oort constants found arc in good agreement with their determinations in other 
works. For example, using data on young objects of the Galactic disk, Bobylev et al. (2008) 
found A = 15.5 ± 0.3 km s"^ kpc"^ and B = -12.2 ± 0.7 km s"^ kpc-^ 

4.1 Parameters of the Velocity Ellipsoid 

We clearly see from Table 2 that the kinematic parameters differ for different groups of giants 
located near the Galactic plane at |Z| < 200 pc (column 4) and at large elevations above the 
Galactic plane \Z\ > 200 pc (column 5). The vertex deviation reaches l^y — 9.1 ± 0°.5 for 
the first group and is a factor of 2 smaller, 4.9 ± 0°.6, for the second group. Note that the 
difference in kinematic parameters is enhanced with increasing elevation above the Galactic 
plane. On the whole, our results are in agreement with those of Rybka (2008). 

From Table 3, we see an increase in velocity dispersions and a decrease in vertex deviation 
with increasing "age" of the sample of stars. Column no. 1 gives the parameters obtained 
for stars with available radial velocities from all three space velocity components. Using the 
radial velocities affects the determination of the parameters of the velocity ellipsoid, with the 
most prominent parameter being li ~ 11°. 7. For other samples, there is also a difference in 
vertex deviation. Thus, for example, h — 5°.l for the sample of oldest stars is considerably 
smaller than li = 9°. 4 for the youngest group. 

On the other hand, since RGC stars with various ages are well mixed in space, using 
various constraints on the spatial coordinates of stars does not lead to any noticeable change 
in the parameters of the velocity ellipsoid (Table 4). 

In contrast to the approach used by Rybka (2008), we traced the changes in the param- 
eters of the velocity ellipsoid with \Z\ (Fig. 2). As expected, our results show that younger 
stars with a lower velocity dispersion are concentrated to the Galactic plane to a greater 
extent. 

4.2 The Contraction Effect 

According to Table 2, the Oort constant K reaches its maximum value for the sample of 
distant (500-1000 pc, d = 0.65 kpc, column no. 4) RGC stars located near the Galactic plane 
(|Z| < 200 pc), —5.3 ±1.3 km s~^ kpc~^, while the contraction velocity is Kd — —3.5 ± 0.9 
km s~^; this value is indicated in Fig. 4 by the filled circle. The open circles in Fig. 4 
indicate five values of Kd from the cited data for which estimates of the average distance are 



available. Below, we list them, as denoted by the numbers in Fig. 4, in order of increasing 
d: 

(1) the contraction velocity determined from the space velocities of A0-A5 giants 
(Bobylev et al. 2006b) for which K = -13.1 ±2.0 km s"^ kpc"^ was obtained, with d = 0.19 
kpc; 

(2) the estimate from Rybka (2004b) obtained using a sample of pallOOO Tycho-2 G5-K0 
giants with |6| < 50°, d > 250 pc, d = 0.364 kpc; 

(3) the estimate from Rybka (2004a) obtained using a sample of ? 9000 main-sequence 
Tycho-2 B0-B9 stars, d = 0.53 kpc; 

(4) the estimate from Rybka (2008) obtained using ^30000 Tycho-2 RGC stars at |6| < 
30° and d > 100 pc for which we found the average photometric distance based on the 
catalog by Gontcharov (2008), d — 0.54 kpc; 

(5) the result obtained using open star clusters with an average age of fslOO Myr (Bobylev 
et al. 2007) at d = 1.05 kpc. 

Note also the result of analyzing OB stars located in the range of distances 600-2000 
pc from Torra et al. (2000) for which K = —2.9 ± 0.6 km s~^ kpc~^ was found. Taking 
d 0.8 kpc as a realistic average distance for these stars (which is not given in Torra et al. 
(2000)), we will obtain Kd — —2.3 km s~^, a value consistent with the data in Fig. 4. A 
comprehensive overview of the results of determining the negative K effect by other authors 
can be found in Fernandez et al. (2001). 

As can be seen from Fig. 4, on the whole, there is good agreement between the various 
contraction velocity estimates. The main purpose of our comparison is to show that the 
results of different authors obtained by different methods agree in sign, i.e., the contraction 
effect actually takes place, and, given the errors, they also agree satisfactorily in magnitude. 
However, the difference for the A0-A5 giants of the OSACA database (Bobylev et al. 2006b) 
is most pronounced. This difference can be explained by the fact that the sample of A0-A5 
stars was specially formed by a kinematic indicator by the maximum K effect. At the same 
time, this result is in good agreement with the data of Table 1, according to which we have 
Kd = -1.8 ±0.5 km s'^ 

What is the maximum contraction velocity? In the opinion of Rybka (2004a), Kd = 
—8.2 ± 0.9 km s~^ at d = 1 kpc. However, this estimate was obtained by extrapolating the 
data on main-sequence dwarfs {d < 0.6 kpc). Analysis of the kinematics of open star clusters 
(Bobylev et al. 2007) showed the magnitude of the velocity Kd to decrease at d > 1 kpc. 

Note also the ambiguity of treating K = 0.5(Mii ±M22) as the effect of contraction in the 
plane. Indeed, as we see from Table 1, M22 = du2/dr2 does not differ significantly from zero, 
with large values of K and Ixy being reached solely through the gradient Mn = dui/dri. 
Thus, we can only assert that noticeable deviations from circular motions are observed. The 
nature of these deviations requires a further study. 

4.3 The Problem of the Galactic Disk Warp 

Rotation eiround the X axis. As can be seen from Table 2 (the third column), — 
— 1.4±0.2 km kpc~^ and M^^ — ±1.1±0.3 km s~^ kpc~^. The negative sign of means 
that the rotation is from the Z axis to the Y axis, i.e., the stars approach to the Galactic 
plane. They approach from the North Pole in the first and second quadrants and from the 
South Pole in the third and fourth quadrants. By analogy with the rotation around the 



Z[Qz = B — A) axis, the angular velocity of rotation around the X axis can be determined 
as Qx = M- - M+ = -2.5 ± 0.3 km s"^ kpc^^ 

This is inconsistent with the analysis of the motions of Ri4000 giants of various spectral 
types 0-M performed by Miyamoto et al. (1993), who found a positive direction of rotation 
around the X axis (directed toward the Galactic center). Note that, in this case, the stellar 
proper motions being analyzed were determined in the FK5 system, which is noticeably 
distorted by the uncertainty in the precession constant. On the other hand, having analyzed 
the proper motions of 0-B5 stars with the Hipparcos (19977) system, Miyamoto and Zhu 
(1998) also reached the conclusion about a positive rotation around the X axis. In our 
notation, Miyamoto and Zhu (1998) obtained fix = SM^^ = 2M+ = +3.8 ±1.1 km s'^ 
kpc^^. However, it should be noted that Miyamoto et al. (1993) and Miyamoto and Zhu 
(1998) used a model where the variables M^^ and could not be separated and, therefore, 
they were assumed to be equivalent. In our approach, which consists in solving Eqs. (7)-(8), 
we can separate these variables. This is possible, because the distribution of RGC stars over 
the celestial sphere is fairly uniform. 

The parameters found with a constraint on |Z| and given in column no. 4 of Table 2 are 
also of great interest. As we see, for the sample of stars close to the Galactic plane (this 
case corresponds largely to the distribution of 0-B5 stars), the solutions for M with index 
3 differ markedly from those obtained for the samples filling uniformly the celestial sphere 
(Table 2, columns nos. 1, 2, 3, 5). 

As can be seen from Table 2 and Fig. 5, the values of the parameters Mg^ and do 
not depend strongly on the distance. However, the motions of nearby stars can be distorted 
by local peculiarities (in particular, by the kinematic peculiarities of the Gould Belt with a 
characteristic radius of about 500 pc). Therefore, the values of these parameters determined 
from distant stars seem more reliable. Thus, we interpret Qx found from relatively distant 
RGC stars as the rotation related to the warp of the Galactic stellar-gaseous disk. 

Note that owing to the use of currently available data that realize the inertial frame of 
reference more accurately than FK5 and the possibility of separating the components 
and when solving Eqs. (7)-(8), wc have every reason to believe that our results are 
more reliable than those of Miyamoto ct al. (1993) and Miyamoto and Zhu (1998). 

Rotation around the Y cixis. There is a large uncertainty with regard to the mag- 
nitude and direction of the rotation around the Y axis related to the disk warp. Thus, 
for example, having analyzed the motions of Ri4000 giants of various spectral types O- 
M, Miyamoto et al. (1993) found a positive direction of rotation around the Y axis: 
fly = 2Mj]^ = — 2M+ = +5.6 ± 1.0 km s^^ kpc~^. Based on an analysis of the Hippar- 
cos data, Drimmel et al. (2000) obtained a completely different estimate in the form of 
precession of OB stars in the ZX plane with an angular velocity of —25 km s~^ kpc~^ 
(i.e., rotation around the Y axis). However, Bobylev (2004) showed that this value could 
be explained by inaccurate referencing of the ICRS/Hipparcos realization to the system of 
extragalactic sources. 

Based on the data of column no. 2 in Table 2, we have fly = M~^—M^^ — —1.5+0.3 km s"-*^ 
kpc~^. As we see, the value of this component is determined almost entirely by M]]^. In this 
respect, M;^/4.74 = —0.41 + 0.05 mas yr"-*^, which is in excellent agreement with the residual 
rotation of the ICRS relative to the inertial frame of reference, = —0.37 + 0.04 mas yr~^, 
estimated by Bobylev and Khovrichev (2006) based on the Tycho-2 and UCAC2 (Zacharias 
2004) Catalogues, is of great interest. Therefore, we cannot say what the contribution related 
to the disk warp is in this case, since these two effects cannot be separated. Obviously, if we 



apply the correction pa —0.4 mas yr^^ (reducing the Hipparcos system to the inertial 
frame of reference), then the "remainder" of the rotation around Y axis will be close to zero. 

On the whole, we may conclude that the kinematic effect related to the disk warp man- 
ifests itself mainly as the rotation around the X axis (see the previous section). The small 
magnitude of the rotation velocity around the Y axis is in agreement with the fact that 
the line of nodes of the HI layer is known to be close to the direction of the Galactic cen- 
ter/anticenter. 

CONCLUSIONS 

We analyzed the proper motions of 95633 RGC stars from the Tycho-2 Catalogue using the 
linear Ogorodnikov-Milne model. The photometric distances to these stars were estimated 
by Gontcharov (2008) based on 2MASS photometric data. 

Based on a subsample of 3632 RGC stars for which not only the proper motions but 
also the radial velocities, photometric distances, and Hipparcos trigonometric parallaxes are 
available, we showed that, apart from the star centroid velocity components relative to the 
Sun (Xq, Yq, Zq), only those parameters of the Ogorodnikov-Milne model that describe the 
stellar motions in the XY plane differ significantly from zero. This allowed us to study the 
plane K effect based the tangential velocities of a considerably larger sample of RGC stars. 

To ascertain the dependence of the kinematic parameters being determined on the dis- 
tances, including the measurement errors in the distances, we performed statistical Monte 
Carlo simulations. Our simulations showed that, first, a significant bias is observed only for 
the velocity 1^0 at large distances and, second, very close {d < 200 — 250 pc) stars should be 
excluded to reliably determine all parameters of the Ogorodnikov-Milne model, except the 
velocity components of the Sun, which, on the contrary, are determined more reliably from 
samples of nearby stars. 

Analysis of the entire sample of RGC stars showed that the Oort constants that describe 
the rotation around the Z axis arc A = 15.9 ±0.2 km s~^ kpc^^ and B = —12.0 ±0.2 km s^^ 
kpc^^. This gives an estimate of the rotation velocity of the solar neighborhood, Vq = 222±2 
km s~^, for Rq — 8.0 kpc and the period of its revolution around the Galactic center, 221 
Myr. 

Judging by the velocity dispersions found, the RGC sample includes stars of various 
ages. We made an attempt to roughly separate the RGC stars by ages based on the parameter 
A^*. It showed that only the fraction of relatively young RGC stars, which accounts for about 
20% of the entire sample, is identified most rehably. 

The most interesting results were obtained using stars located in the range of distances 
500-1000 pc. The RGC stars located near the Galactic plane were shown to have a noticeable 
contraction effect (observed in the XY plane) that depends on the heliocentric distances of 
the stars. Thus, for example, for the sample of RGC stars located near the Galactic plane 
(|Z| < 200 pc), the contraction velocity is Kd — —3.5 ± 0.9 km s~^ and the vertex deviation 
is Ixy — 9.1 ± 0°.5. For stars located well above the Galactic plane (|Z| > 200 pc), these 
effects are less pronounced, being Kd = —1.7 ± 0.5 km s~^ and l^y = 4.9 ± 0°.6. 

Based on RGC stars with distances in the range 500-1000 pc, we found the angular 
velocity of rotation around the Galactic X axis to be —2.5 ± 0.3 km s^^ kpc~^, which 
we interpret as the rotation related to the warp of the Galactic stellar-gaseous disk. The 
direction of this rotation shows that the stars approach the Galactic plane. 
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Table 1: Kinematic parameters obtained by solving the system of equations (3)-(5) 







A 












2 






km s^ 


1 


9.9 ±0.4 


9.5 ±0.4 






km s~ 


1 


16.8 ±0.4 


16.3 ±0.4 






km s~ 


1 


7.3 ± 0.4 


6.9 ±0.4 


Mil, 


km 


^ kpc" 


1 


-15.5 ±3.3 


-14.3 ±3.7 


Mi2, 


km s" 


^ kpc" 


1 


28.7 ±3.2 


29.3 ±3.5 


Mi3, 


km s~ 


^ kpc" 


1 


-2.8 ±2.6 


0.6 ±2.9 


M21, 


km s~ 


^ kpc" 


-1 


8.6 ±3.3 


5.9 ±3.7 


M22, 


km s~ 


^ kpc" 


-1 


-1.3 ±3.2 


-2.2 ±3.5 


M23, 


km s~ 


^ kpc" 


1 


-5.8 ±2.6 


-3.8 ±2.9 


M31, 


km s~ 


^ kpc" 


1 


-6.7 ±3.3 


-4.3 ± 3.7 


M32, 


km s~ 


^ kpc" 


-1 


-5.3 ±3.2 


-4.1 ±3.5 


M33, 


km s~ 


^ kpc" 


-1 


-1.1 ± 2.6 


-4.3 ±2.9 


^, 


km s~ 


^ kpc" 


1 


18.6 ±2.3 


17.6 ±2.5 




km s~ 


^ kpc" 


1 


-10.0 ±2.3 


-11.7 ±2.5 


c, 


km s~ 


^ kpc" 


-1 


-7.1 ±2.3 


-6.1 ±2.5 




km s~ 


^ kpc" 


1 


-8.4 ± 2.3 


-8.2 ±2.5 






d, kpc 


0.22 


0.22 



Note. The parameters presented in columns nos. 1 and 2 were obtained using trigonometric 
parallaxes and photometric distances, respectively; N^, is the number of sample stars. 



Table 2: Kinematic parameters obtained by solving the system of equations (7)-(8) 







d < 1000 


PC 


600 


- 1000 


PC 


700 


- 1000 pc 


600 


- 1000 pc 


500 - 1000 pc 
























\Z\ < 200 pc 


|Z| > 200 pc 








95633 




50595 




19200 




22684 


27911 




1 


2 


3 


4 


5 




7) , km s ^ 


7 


84 


± 


09 


8 


55 ± 


13 


8 


67 ± 0.22 


8 


49 ± 0.18 


8.52 ± 0.18 


Y 


2 , km 


15 


94 


± 


10 


17 


74 ± 


14 


19 


19 ± 0.26 


16 


14 ± 0.23 


18.48 ± 0.19 


z 


7) , km s~ ^ 


6 


60 


± 


08 


7 


17 ± 


12 


7 


59 ± 0.21 


6 


86 ± 0.15 


7.59 ± 0.19 


M+(A), 


km/s/kpc 


15 


87 


± 


20 


15 


86 ±0 


24 


15 


61 ±0.37 


16 


25 ± 0.35 


15.33 ±0.35 




km/s/kpc 


-1 


40 


± 


17 


-1 


41 ± 


21 


-1 


80 ± 0.31 


1 


83 ± 0.95 


-1.70 ± 0.26 




km/s/kpc 


-2 


01 


± 


18 


-1 


93 ± 


22 


-1 


97 ± 0.32 


-6 


28 ± 0.77 


-1.70 ± 0.27 


M (B), 
<i-22. 


km/s/kpc 


-11 


99 


± 


15 


-11 


99 ± 


18 


-12 


18 ± 0.26 


-11 


82 ± 0.24 


-12.26 ± 0.27 


km/s/kpc 


-7 


86 


± 


38 


-7 


43 ± 


45 


-6 


47 ± 0.68 


-10 


73 ± 0.62 


-5.29 ± 0.66 


M+ 


km/s/kpc 


-0 


45 


± 


23 


-0 


47 ± 


27 





15 ± 0.40 


-3 


90 ± 0.82 


-0.31 ± 0.38 




km/s/kpc 





98 


± 


22 


1 


13 ± 


26 


1 


49 ± 0.38 


-1 


79 ± 1.01 


0.88 ± 0.36 


M 

33-22' 


km/s/kpc 


-1 


36 


± 


43 


-0 


99 ± 


51 


-0 


47 ± 0.75 


-0 


02 ± 1.27 


-0.27 ± 0.61 




d, kpc 





52 


± 


21 





68 ± 


14 





81 ± 0.11 





65 ± 0.12 


0.70 ± 0.15 


C, km 


~^ kpc^l 




-3.9 ± 0.2 




-3.7 ± 0.2 




-3.2 ± 0.3 




-5.4 ± 0.3 


-2.7 ± 0.3 


K, km i 


-1 kpc^l 




-2.6 ± 0.5 




-2.7 ± 0.6 




-2.8 ± 0.8 




-5.3 ± 1.3 


-2.4 ± 0.7 


K ■ 


d, km s ^ 




-1.4 ± 0.3 




-1.8 ± 0.3 




-2.3 ± 0.7 




-3.5 ± 0.9 


-1.7 ± 0.5 








7.0 ± 0.3 




6.6 ± 0.4 




6.9 ± 0.6 




9.1 ± 0.6 


4.9 ± 0.6 



Table 3: Parameters of the residual velocity ellipsoid 



AT 


ODOZ 


oyy (y 


oDOzy 




s 




1 


2 


3 


4 


5 


0-1 


32.1 ±0.4 


29.1 


35.6 


19.7 


0.50 


0-2 


21.2 ±0.4 


21.7 


25.9 


13.1 


0.40 




17.9 ±0.3 


16.8 


20.6 


10.8 


0.25 


^1 


11°.7±4°.0 


6°.3 


5°.l 


9°.4 


0°.4 


h 


101°.7± 1°.4 


96°. 3 


91°.8 


100°. 1 


0°.4 


h 


287°.6± 1°.4 


260° 


253° 


265° 


15° 


hi 


-4°.l ±0°.8 


1°.0 


0°.8 


3°.0 


0°.2 




3°.9±0°.8 


3°.6 


5°.0 


11°.9 


0°.6 


h 


86°. 1 ± 1°.8 


86°. 3 


87°. 9 


77°. 7 


0°.5 



Note. Columns nos. 1 and 2 4 give the parameters determined from the stellar space velocities 
and only from the stellar proper motions, respectively; column no. 2 for < AS, column no. 3 for 
—0.38 < AS < 0, column no. 4 for —1 < AS < —0.38; column no. 5 gives a Monte Carlo estimate 
of the error e. 



Table 4: Parameters of the residual velocity ellipsoid obtained using the proper motions of RGC 
stars for the samples of Table 2 





d < 1000 pc 


500 - 1000 pc 


700 - 1000 pc 


500 - 1000 pc 


500 - 1000 pc 










\Z\ < 200 pc 


\Z\ > 200 pc 




95633 


50595 


19200 


22684 


27911 




1 


2 


3 


4 


5 


0-1 


30.1 


32.2 


35.0 


31.4 


33.0 


0-2 


21.7 


23.2 


25.3 


20.3 


24.2 




17.2 


18.7 


20.7 


17.3 


20.4 


h 


6°.0 


5°.2 


5°.5 


6°.6 


3°. 8 


h 


96°.0 


95°. 2 


95°.5 


96°.6 


93°.8 


h 


257°. 9 


91°. 1 


83°. 7 


242°. 7 


80°. 1 


hi 


0°.7 


-3°.0 


-1°.2 


0°.5 


-l°.l 


h2 


2°.3 


-4°.2 


-5°.7 


0°.7 


-4°. 6 


hs 


87°. 6 


85°.8 


84°. 2 


89°. 1 


85°. 3 




Fig. 1. Distributions of the samples of RGC stars in the XY, XZ, and YZ planes (the left, 
central, and right columns, respectively): the complete sample (a), the stars with distances 
d < 400 pc (b), 400 pc< d < 600 pc (c), and 600 pc< d < 1600 pc (d). 
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Fig. 2. Principal semiaxes of the residual velocity ellipsoid a^, a2, versus \Z 




Fig. 3. Smoothed distributions of UV velocities: (a) sample (1) for < AS; (b) sample (2) 
for -0.38 < AS <Q; (c) sample (3) for -1 < AS < -0.38 (see the text). 
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Fig. 4. Velocity Kd versus heliocentric distance d\ the open circles mark the cited data (see 
the text) and the filled circles indicate the value of Kd from Table 2. 




Fig. 5. Parameters of the Ogorodnikov-Milne model versus distance calculated without 
introducing any model errors. 




Fig. 6. Parameters of the Ogorodnikov-Milne model calculated with model distances versus 
distance; the random errors in the individual distance of a star are 20%. 




Fig. 7. Parameters of the Ogorodnikov-Milne model calculated with model distances versus 
distance; the random errors in the individual distance of a star are 60%. 



